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Abstract 

The exclusive production cross sections for iv and (p mesons have been mea- 
sured in proton-proton reactions at piab = 3.67 GeV/c. The observed (j)/uj 
cross section ratio is (3.8 ± 0.2lQ;g) x 10^^. After phase space corrections, 
this ratio is enhanced by about an order of magnitude relative to naive pre- 
dictions based upon the Okubo-Zweig-Iizuka (OZI) rule, in comparison to an 
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enhancement by a factor ~ 3 previously observed at higher beam momenta. 
The modest increase of this enhancement near the production threshold is 
compared to the much larger increase of the ^/w ratio observed in specific 
channels of pp annihilation experiments. Furthermore, differential cross sec- 
tion results are also presented which indicate that although the (f) meson is 
predominantly produced from a ^Pi proton-proton entrance channel, other 
partial waves contribute significantly to the production mechanism at this 
beam momentum. 
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Typeset using REVT^ 



I. INTRODUCTION 



The proton is a complex composite object in which sea-quarks and gluons make signif- 
icant contribution to the structure functions. Investigation of the strange structure of the 
nucleon is currently of great interest since strange quarks appear only as sea quarks in most 
models of the nucleon in contrast to the light quarks {u and d) that appear as both valence 
and sea quarks. 

An important experimental approach to study the role of strange quarks in the proton's 
wave function is to measure the relative production of and u vector mesons. Since the 
singlet-octet mixing angle of the ground state vector meson nonet is close to the ideal value 
(i.e. tanOideai = 05 the u meson consists almost completely of light valence quarks 

and the (p meson almost completely of strange valence quarks. According to the Okubo- 
Zweig-Iizuka (OZI) rule, processes with disconnected quark lines in the initial or final state 
are suppressed [0-^]. As a result, the production of the meson is expected to be strongly 
suppressed relative to the u meson in hadronic reactions with no strange quarks in the initial 
state. Under the assumption that the OZI rule is exactly fulfilled, the meson can only be 
produced by the small admixture of light quarks to its wave function. Thus, the relative 
production cross sections for the and u mesons (-R^/,^) can be calculated by the following 
formula: 

R^/u; = / X tan^ 5 = / X 4.2 ■ 10~^ (1) 

where / is a correction for the available phase space, and 5 ~ 3.7° is the deviation from the 
ideal mixing angle 0. 

In certain hadronic reactions it has been observed that the ratio of the exclusive cross 
sections for the and u meson production reactions significantly exceeds estimates based on 
simple quark models (see e.g. 0). For pp annihilation at rest this enhancement is observed to 
be particularly dramatic ||8HlO|. This apparent violation of the OZI rule has been interpreted 



as evidence for a non-negligible negatively polarized ss Fock component to the proton's wave 
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function 0,0. Here it is important to distinguish between extrinsic and intrinsic quark and 



gluon contributions to the nucleon sea [0,0 , since other explanations stress the importance 
of higher order rescattering processes |T^-|19|, thereby "avoiding" the OZI rule. 

Recently, there have been numerous experimental investigations on the contribution of 
strange quarks to the proton: At low momentum transfers there are measurements of 
parity violating polarized electron-proton scattering to determine the strange electromag- 
netic form factors of the nucleon. First results indicate very small effects [pm plj, and thus 
the scale of the strange vector current (A^|s7^s|A^) can not yet be determined (see e.g. [0). 
Spin dependent structure functions have also been measured in a number of inclusive [^^-|3l|| 



and semi-inclusive p^j33[] experiments on deep inelastic scattering of polarized leptons from 



nucleons, suggesting that the strange sea quarks are weakly polarized against the proton 
helicity. Furthermore, it has commonly been discussed that the magnitude of the St^at term 
in low energy pion nucleon elastic scattering indicates a possibly sizable, but very 

uncertain scalar density of strange sea quarks (A^|ss|A^) in the nucleon's wave function. 

The negatively polarized intrinsic strangeness model was in part motivated by the strong 
correlation of the meson yield to the spin triplet fraction of the annihilating pp system . 



Since meson production in pp reactions at threshold must proceed via the spin triplet 
entrance channel, further insight on the origin of the enhanced production in pp could be 
provided by studying near-threshold pp reactions, where predictions based upon the strange 
sea quarks in the nucleon and meson rescattering models might be expected to differ. 

To address this problem we have measured the production of and uj mesons in proton- 
proton reactions with the DISTO spectrometer ||3^ at the Saturne accelerator in Saclay. We 



present here total and differential cross sections that were determined for a beam momentum 
of 3.67 GeV/c which corresponds to a much lower available energy above the production 
threshold (Q = ^ - ^ = 83 MeV) than the existing data (Q > 1.6 GeV) |0|-|3l. In 



previous letters P^ , ^ we have presented the first results close to the meson production 
threshold in proton-proton reactions, where it was determined that the (p/uj cross section 
ratio increases slightly toward threshold. Since these publications we have made several 



improvements to the data analysis, including refinements of the acceptance correction, im- 
provements in the tracking algorithm and calibrations, as well as the analysis of a factor four 
more statistics. These results are presented here in more detail, together with differential 
cross section distributions that were previously not available. All of the previously published 
data are consistent with the results presented here. 

This paper is organized as follows: in the next section experimental details of the detector 
and data analysis are presented. The resulting and uo meson total and differential cross 
sections are presented in the ensuing section, followed by a discussion of these results and a 
summary. 



II. EXPERIMENT 



A. Apparatus 



A proton beam from the SATURNE synchrotron with momentum piab = 3.67 GeV/c 
was directed onto a liquid hydrogen target of 2 cm length, and multiple charge particle final 



states were measured with the DISTO spectrometer |39|, which is shown schematically in 
Figure 0. Charged particles were tracked through a magnetic spectrometer and detected by 
a scintillator hodoscope and an array of water Cerenkov detectors. The magnetic spectro- 
meter consisted of a dipole magnet (1.0 T-m), 2 sets of scintillating fiber hodoscopes inside 
the field, and 2 sets of multi wire proportional chambers (MWPC) outside the field. The 
large acceptance of the spectrometer ±15° vertical, ±48° horizontal) allowed for coinci- 
dent detection of four charged particles, which was essential for the kinematically complete 
reconstruction of many final states (ppTr+vr^, pp7r+7r^7r°, ppK~^K~, pKA, pKH). Particle 
identification and 4-momentum conservation served as powerful tools for background rejec- 
tion. Event readout was triggered by a multiplicity condition on the scintillating fiber and 
hodoscope detectors, selecting events with at least three charged particles. Details of the 



experimental apparatus can be found in Ref. |39 . 
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B. Data Analysis 



The exclusive reactions pp pprj and pp ppoo were identified via the vr'^Tr^Tr^ decay of 
the T] and lj mesons with partial widths T/Ttot = 0.232 and T/Ttot = 0.888, respectively [|l|. 
For the pp ppcf) reaction, the was observed via its K^K~ decay with T/Ttot = 0.491. 
Since these channels all have four charged particles in the final state (three positive and 
one negative), these events could be reconstructed by applying particle identification and 
kinematical constraints on an event-by-event basis to the same sample of four track events. 

1. Tracking 

After calibration, the data from the magnetic spectrometer consist of a collection of 
points in space determined by the various detector components. Thus, the first major 
step in the data analysis is to determine the trajectories of the charged particles from the 
detector position information. In general, a particle's trajectory is defined by 5 independent 
parameters. These parameters can be chosen in numerous ways, however they must span 
the space of allowed possibilities. The choice of track parameters used in this experiment to 
define the trajectories was motivated by the dipole magnetic field geometry, and was based 
on a right handed coordinate system with its origin in the center of the target, z in the 
beam direction at the target, and y (vertical) along the magnetic field direction. The five 
track parameters used were X and Y to describe the intersection of the track with the x — y 
plane through the origin, to describe the angle of the track's projection onto the x — z 
plane at the point {X, Y,0), m to describe the ratio of the particle's momentum along y to 
the momentum in the x — z plane, and Pxz for the reciprocal of the momentum in the x — z 
plane. 

The observed positions in the detectors can be calculated as a function of these five 
parameters by detailed Monte-Carlo simulations of the detector performance. The goal of 
the tracking algorithm is to invert this function, i.e. determine the track parameters from 
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the measured position information. This is done in a two step procedure: In the first step 
(track-search) it is determined which hits belong together to form a track. Then, in the 
track-fit procedure the five track parameters are determined from the hits in an iterative, 
quasi-Newton procedure by interpolating within a matrix of reference trajectories. 

2. Particle Identification 

Particle identification was achieved using the light output from the water Cerenkov 
detectors, which provided good 7r+ - proton separation spanning a wide range of mo- 
menta, as well as identification in a restricted range above the kaon Cerenkov thresh- 
old j)K,th = Mxc/y/n"^ — 1. Water was chosen as the Cerenkov radiator {n = 1.33, 
PK,th = 560 MeV/c) in order to match the momentum range of the kaons from the reaction 
pp ppK^K~ , which are distributed around 700 MeV/c. 

The inclusive correlation of Cerenkov amplitude versus momentum is shown in Figures 0a 
and ^ for negatively and positively curved tracks, respectively. These figures have a loga- 
rithmic intensity scale with a factor 2.4 between intensity contours. In Figure ^ there is a 
clear band from the vr" mesons. The K~ mesons are not visible in this inclusive distribution 
because of the ~ 10^ larger pion yield. On the other hand, there is already a weak indication 
of the mesons in Figure ^jc in addition to the prominent vr^ meson and proton lines. The 
relative yield of a given kaon species can be enhanced by imposing two requirements on the 
data set: The first is that the oppositely charged meson is consistent with being a kaon, and 
the second is 4-momentum conservation for the given event hypothesis (i.e. pp — >■ ppK^K^). 
The remaining yield after imposing these restrictions is shown in Figures 03 and where 
signals from the and K~ mesons are now clearly visible. 

3. Event Reconstruction 

The data sample analyzed here with four charged particles contained events from many 
different types of reactions. The individual reactions were identified by a combination of par- 
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tide identification of tlie cliarged particles in tlie final state and kinematic constraints on the 
measured particles. The charged particles in the final state, predominantly protons, pions, 
and kaons, were identified with the water Cerenkov detectors as described above. How- 
ever, since the same final state could be produced via several reaction channels, kinematic 
conditions were required to discriminate among the particular reaction channels. 

All reaction channels studied here contained at most one unmeasured particle in the 
final state, thus the measured four track events were kinematically complete. As a result 
4-momentum conservation was imposed as a powerful constraint to select events of a given 
hypothesis. For this the invariant mass and missing mass of various particle configurations 
were calculated: The invariant mass (Mj„i,) is the total energy in the reference frame of a 
given number of observed particles. Mj„t, can be calculated from the total energies {E) and 
momenta (p) of the n individual particles using the following formula: 



The missing mass analysis was used to determine the rest mass of an unobserved recoil 
particle. The missing mass {Mmiss) is given by the following formula: 



where the index i runs over all m observed particles. The specific application of these 
constraints is discussed in the following sections. 

The results presented here on the production of the meson were determined by selecting 
events of the type pp — » ppK^K~ . For this event hypothesis, the proton-proton missing 
mass (M^^g) must equal the invariant mass of the kaon pair (M^^). Since all four particles 
in the final state have been measured, the event reconstruction has been performed requiring 
that the four particle missing mass [M^^^) be equal to zero. The distribution of {M^^f- 
(M^igg)^ is plotted in Figure ^ including kaon identification (solid histogram) based on the 
Cerenkov detectors. The peak at {M^^Y — (M^^^^)^ ^ results from events consistent 
with the pp ppK^K~ hypothesis, and is superimposed on a background resulting from 




(2) 




(3) 
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imperfect tt — K separation in the Cerenkov detectors in a small fraction of events of the 
type pp — s> pK^A ppK^ix^ or pp —* ppn^Ti^X . The dashed histogram is an estimate of 
the background by scaling the inclusive distribution without kaon Cerenkov requirements by 
a factor 0.002 in order to match the data in Figure ^ above 0.15 GeV^/c^. The background 
from non ppK^K~ final states comprises only about 2.2% of the yield in the range | [M^^Y — 
{M^^gg)^\ < 0.09 GeV^/c^, which marks the range where events were accepted for the further 
analysis. 

The exclusive reaction channel pp — > ppuj was measured via the uj —>■ Tr+Tr^Tr" decay 
and was selected by a missing mass analysis since the 7r° decays primarily into two pho- 
tons (r^^/r^Q^ = 0.988) which are not detected in the spectrometer. Events of the type 
pp-K^-K^-K^ were selected by first requiring that the Cerenkov amplitudes associated with the 
four charged tracks be consistent with the hypothesis ppn^TT^. Furthermore, background 
from four body reactions of the type ppn^n^ with M^^^ ~ is partially suppressed by 
requiring the invariant mass of the pion pair to fulfill {M[^^'^ )^ < (M^^ — M^o)^. After 
these requirements, the distribution of the four particle missing mass squared (M^^^^ )^ 
versus the proton-proton missing mass (M^^^^)^ is shown in Figure ^ with a linear scale 
for the intensity contours. As seen in this figure, there is a very strong signal for the 
pp ppuj ppn+TT-n° reaction at {M^^^'^^'f ^ (M^o)^ and (MZss)^ ~ (M^)^. The 
projection of this distribution onto the (M^^^^)^ axis is shown in the lower frame of Figure ^ 
with the additional requirement that (M^^j^'^ )^ ^ (M^o)^. In this figure, there are clear 
peaks from the tt+tt^tt*^ decay of the rj and u mesons. The structure at very low (M^^^)^ is 
due to contamination of four body events of the type ppn'^TT' . 



C. Acceptance Correction 

The relative acceptance of the apparatus for the ppt], ppuj and ppK~^ production chan- 
nels has been evaluated by means of Monte Carlo simulations, which after digitization, were 
processed through the same analysis chain as the measured data. The relative acceptance 
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could be determined independent of the actual phase space distribution of the particles in 
the final state, because the following two requirements are fulfilled with the DISTO spec- 
trometer: (i) the detector acceptance was determined as a function of all relevant degrees 
of freedom in the final state, and (ii) after accounting for the azimuthal and refiection sym- 
metries, the detector acceptance was non-zero over the full kinematically allowed region. 

In order to determine the detector acceptance, the relevant degrees of freedom were each 
subdivided into a specific number of bins, thus defining a multi-dimensional matrix. Next, 
the number of both the generated and the reconstructed events from the simulations were 
stored in separate copies of this matrix. Finally, the bin-by-bin ratio of the number of 
generated to the number of reconstructed events provided the efficiency correction matrix. 

When analyzing the data the acceptance correction procedure was performed on an 
event-by-event basis by first determining the phase space bin of the particular event and 
then incrementing the distributions by the corresponding acceptance correction factor from 
the acceptance correction matrix described below. These correction factors include the 
appropriate partial widths and the simulations account for the lifetime and different decay 
modes of all unstable particles involved. 

1. ppK^K~ final state 

Although there are in general 16 degrees of freedom for 4 particles in the final state, 
numerous constraints exist which significantly lower the total number of independent degrees 
of freedom. For instance, particle identification and 4-momentum conservation reduce the 
number of degrees of freedom to 8 for the ppK^K~ final state. The remaining degrees of 
freedom can be parameterized in numerous ways, with the only restriction that the variables 
chosen must span the space of allowed possibilities. 

Without loss of generality the eight kinematic degrees of freedom for the ppK~^K~ final 
state can be parameterized as if the reaction proceeded in two steps, i.e. pp — > P1P2X — > 
PiP2K'^K~. In this case, five variables can be used to describe the three body system ppX, 



10 



and the remaining three degrees of freedom are then required to uniquely determine the 
decay of the intermediate state X. This choice of kinematic variables provides a complete 
basis and does not require that the reaction actually proceeds via the intermediate state 
X. The ppX system is parameterized by the two invariant mass combinations (Mf^^)"^ and 
{Mf^^Y (i.e. Dalitz plot variables of the ppX system), and three Euler angles to describe 
the orientation of the ppX decay plane in the center of mass reference frame (CM): the 
polar and azimuthal angles of the X intermediate state Qcmi 'Pcmi ^ rotation of the 
ppX decay plane around the direction of X, iI^cjm- Furthermore, three variables are required 
to describe the decay of X: the mass Mx = , and the emission angles of one of the 
kaons in the X reference frame and (px where these angles are measured relative to the 
beam direction^. The beam direction is chosen as the reference direction for the X decay 



to facilitate comparison in Sec. piB| with simple expectations for 0-meson decay when it is 
produced near threshold in pp reactions. 

The differential cross section for X production must be isotropic as a function of the az- 
imuthal angle (pcM- Therefore, in the simulations (pcM been integrated with an isotropic 
distribution to determine the detector acceptance as a function of the remaining seven de- 
grees of freedom. 

Due to the high dimensionality of this problem, even a small number of bins per kinematic 
variable result in an acceptance correction matrix that is too large to calculate with the 
available computational facilities. For a subdivision of each degree of freedom into = 10 
bins, and an average of M = 20000 events simulated per bin of the acceptance matrix, a total 



The following conventions are used to label the various angles that appear in this report: O 
corresponds to polar angles measured with respect to the beam axis, and ^ corresponds to polar 
angles measured with respect to any other quantization axis, which will be explicitly stated. Az- 
imuthal angles are labeled by (p and "0! respectively. The subscript denotes the reference frame in 
which the angle in measured, and the superscript denotes the particle being measured. 
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of M X A^^ 2 X 10^^ would have to be simulated. Thus, in order to explore the dependence 
of the detector acceptance on the remaining seven kinematic variables, two subsets of the 
general acceptance correction matrix (Matrix A and Matrix B) have been generated by 
integrating over several of the kinematic variables as discussed below and summarized in 
Table | 

The simulations for Matrix A contain two simplifications in order to make the calculations 
tractable. The first simplification is the requirement that the mass of the intermediate state 
X be equal to the meson mass. Of course this requirement restricts the use of Matrix A 
to the pp(f) reaction only; however, little systematic error is introduced because the shape of 
the mass distribution of the meson is well known and the detector acceptance varies little 
over it. Furthermore, an isotropic angular distribution of the differential cross-section for 
meson production versus Qcm has been assumed for Matrix A. This assumption is consistent 
with the observation that this angular distribution turns out to be isotropic when Matrix 
B (see below), which includes explicitly the Qqm dependence, is used for the acceptance 
correction. Therefore, in the simulations to calculate Matrix A, the intermediate state X 
was given the meson mass and the angular variables Qcm (pQj^j were integrated over 
with isotropic distributions when determining the acceptance as a function of the remaining 
five variables. 

In order to overcome the restrictions on the usage of the acceptance correction Matrix A 
mentioned above, the available degrees of freedom in the ppK^K~ final state were subdi- 
vided differently, as listed in Table |, resulting in the acceptance correction Matrix B. Here 
the variables and were explicitly included in Matrix B. However, an isotropic 

distribution for the decay X — > K^K^ was used in order to reduce the dimensionality of 
the matrix to a solvable level. The validity of using an isotropic decay of X ^ K^K~ can 
be judged by the data presented in Sec. |111B| using Matrix A. Furthermore, the variables 
(^Mfjl^^ and (^M^^^^ were integrated over with a weighting according to three body phase 
space. This was motivated by the observation that the physical distributions show only 
small deviations from isotropy and double checked to give consistent results with a different 
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generator for the ppX system. 



2. pp-K'^ir TT^ final state 

For the measurements of the 1] and u meson production there are 5 particles in the final 
state of each channel observed. The 20 degrees of freedom associated with 5 particles are 
reduced to 12 because of 4-momentum conservation and four particles have been identified. 
Two further degrees of freedom are eliminated by requiring M^^^^ = M^^^, respectively, and 
M^^,^/" = M^o. Of the remaining 10 degrees of freedom, five are related to the ppX system 
where X = ri,uj and five to the corresponding decay X tt+tt^tt^. 

The ppX system is parameterized by the Dalitz plot variables {Mf^^Y and {Mf^^Y, as 
well as three Euler angles to describe the orientation of the ppX decay plane: 6^^, (pcM-i 
and Tp'cjM- discussed above, the differential cross section is isotropic in (p-^Mi ^^us there 
are in effect 4 independent degrees of freedom plus those associated with the decay of the 
meson X. The matrix element associated with the uj tt+tt^tt'^ decay was taken from |^ 



and verified to be consistent with the data from ^7] , ^ , thus allowing the corresponding 
variables to be integrated over. The matrix element for the rj decay was taken from |49H5l|. 
Furthermore, we have assumed an isotropic orientation of the uo decay plane, which was 
verified to be consistent with the data. Thus, four dimensional efficiency matrices were 
calculated for the 7] and uj production reactions, as summarized in Table |I[ 

Although all bins that are kinematically allowed have a finite acceptance, there are some 
phase space bins with very low acceptance. These bins are associated with backward emission 
of the rj or UJ meson (in the CM frame). Because the initial system involves two identical 
particles, the physical distribution must have a symmetry about Qcm = 90°- Thus, in 
order to reduce the systematic error associated with the very large acceptance corrections at 
backward angles, all integrated cross section results for the rj and uj meson production were 
calculated for the forward hemisphere only (and then multiplied by two). The differential 
cross sections however could be extended to backward angles, up to the limits of where the 
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acceptance correction method remained valid, thereby providing an additional check of the 
calculations. 



D. Absolute Normalization 

The absolute cross section normalization was determined by measuring the yield of a 
given channel relative to that of a simultaneously measured channel with known cross section. 
This method to determine the absolute normalization was chosen because it reduced the large 
systematic uncertainty associated with the absolute calibrations of both beam intensity and 
trigger efficiency. For this work the reference channel was the reaction pp — > pprj, for which 
a large amount of existing data [^^-|63| are summarized in the upper frame of Figure ^. 



In order to provide the absolute cross section calibration, the existing published data 
were interpolated to estimate the 1] production cross-section at the beam momentum of the 
present measurement. The solid and dashed curves presented in the upper frame of Figure |^ 
correspond to two different parameterizations of the measured cross section values. The 
solid curve corresponds to a polynomial of sixth order and the dashed curve corresponds to 
a parameterization of the following form: 



(Tpp^pp„ = a (^— - 1 j (^— j • (4) 

where ^/sq = (2Mp+M^) is the CM energy at the t] meson production threshold, and a, b and 
c are free parameters. The vertical dotted line marks the available energy {Q = 0.554 GeV) 
of this measurement. 

Both parameterizations describe the existing t] total cross section data well with the 



exception of the measurement at pbeam = 2.8 GeV/c by E. Pickup et al. ||56|, which is under- 
estimated. In the upper frame of Figure ^ both values cited in Ref . |Q , for identification 
via the 1] —>■ tt+tt^tt'^ and r] neutrals decay channels are plotted at the corresponding 
CM energy above threshold Q = 0.275 GeV. This discrepancy has been neglected since that 
measurement |^ is subject to a large systematic error associated with a quite substantial 
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background subtraction. The average of these two interpolations at Q = 0.554 GeV is our 
estimate of 135 ± 35 fih for the total cross section of the reaction pp — > ppr] at our beam 
momentum. The systematic error from the absolute normalization (±26%) is determined 
from the range of the parameterizations and is similar to or smaller than the systematic 
error from the combination of all other sources (i.e. 32% for the (p/rj ratio, see Table y). 
Furthermore, our estimate is in good agreement with one boson exchange model calculations 



by Vetter et al. 64 , who predict ~ 120 fih at this energy. 



E. Systematic Errors 

Due to the large amount of data collected, the statistical errors are relatively small, 
and the experimental error is dominated by systematic uncertainty and systematic bias. 
This section summarizes the effects studied to estimate the magnitude of the systematic 
uncertainty, as well as systematic deficiencies of the acceptance correction method and trigger 
bias. 

Systematic uncertainty of the results quoted here have been studied in detail for the 
following effects: 

• For the acceptance corrections, the statistical uncertainty of the simulations has been 
studied by comparing the reconstructed particle ratios using acceptance correction ma- 
trices based on different subsets of the simulated data. Furthermore, effects due to the 
finite binning of the kinematic variables, as well as the finite detector resolution have 
been considered. Moreover, an additional cross check of the acceptance corrections 
was provided by the observation that the total (j) meson yield is identical (within the 
error range) when determined using both acceptance correction Matrices A and B. 

• A significant source of uncertainty is related to the meson yield determination (back- 
ground subtraction and peak shape parameterization). This effect was studied by 
varying the parameterization of the line shapes and background, as well as by varying 
the fit ranges. 
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• The four track trigger used to record the data presented here had a shghtly different 
efficiency for ppTT^ir" and ppK^K^ events. Although this effect has been corrected 
for, as discussed below, effects such as noise and cross-talk in the multi-anode photo- 
multipliers used for readout of the scintillating fiber detectors introduce an uncertainty 
to the magnitude of this correction, which we denote as "trigger bias" . 

• Similarly, the tracking efficiency of the data analysis procedure may vary slightly 
between measured and simulated data, primarily due to uncertainty of the actual wire 
chamber efficiency. In particular, the single track efficiency and correlated efficiency 
losses due to small spatial separation between tracks have been studied in detail. 

• Since the data presented here were collected over an extended period of time, residual 
efficiency losses due to long term drifts in the calibration of the electronics may remain. 
This effect was examined by comparing the reconstructed particle ratios in temporally 
separated subsets of the data. 

• Finally, uncertainty arising from the particle identification in the Cerenkov detectors 
and the kinematic conditions has been estimated by comparing the efficiency loss due 
to each individual restriction between the simulated and the measured data. 

These results are summarized in Table |Tl|. The total systematic uncertainty quoted is deter- 
mined by quadratically adding the individual terms applicable for a given measured cross 
section ratio. 

In addition to the systematic uncertainty, several effects lead to a systematic bias of the 
calculated cross section ratios. To estimate the magnitude of the systematic bias we have 
examined the following effects. 

• The liquid hydrogen target was contained in a vessel in which background reactions 
were produced. These "non-target" events could be well separated by the vertex loca- 
tion in those events with complete scintillating fiber information. Since this separation 
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was not possible in the remaining ^ 50% of the events, the relative contamination de- 
termined in those events with full fiber information was subtracted from the complete 
data sample. 

The data collection trigger required that three of the four charged particles must be 
observed in the scintillating fiber detectors. Since the fiber detectors had a different 
response for minimum ionizing pions compared to kaons, which were slower, the trigger 
efficiency was higher for ppK^K~ events than for ppn^Tr— events. The variation of the 
fiber efficiency has been studied as a function of particle velocity. The magnitude of the 
trigger bias was thus determined by combining the average kaon and pion efficiencies 
with the "three of four" trigger condition and an estimate of the cross talk in the fiber 
detectors. 

Finally, in the data analysis, conditions were placed on the Cerenkov amplitude of 
the individual tracks, as well as on several kinematical quantities. The efficiency with 
which a particular event type is accepted by all these conditions is strongly related to 
the detector resolution, which may be imperfectly modeled in the simulations, leading 
to a bias in the efficiency correction matrices. The magnitude of this effect was exam- 
ined for each meson production reaction by comparing the acceptance loss from each 
individual selection criterion in the data to the simulations. 



The individual contributions to the systematic bias are summarized in Table |ITI| as the factor 
by which they modify the given particle cross section ratio. Based on these considerations, 
the reconstructed (p/uj, <f)/i], and 00/ rj cross section ratios have been multiplied by the factors 
0.97, 0.98, and 1.02, respectively, to account for the systematic bias. The correction factor 
for the total K'^K^/rj cross section ratio is the same as for the (p/r] ratio. 
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III. RESULTS 



After full acceptance corrections, the relative yields for exclusive ri,uj, and meson 
production in proton-proton reactions at 3.67 GeV/c have been determined. These results 
are summarized in Table 0, including the statistical and systematic error, and have been 
combined with the absolute rj meson production cross section (see Sect. [11L)| ) to determine 
the absolute cross section results presented below. 



A. io Meson Production 

Using the r] meson yield as the absolute normalization, as discussed above, the total cross 
section for the reaction pp —>■ ppu is determined to be (50 ± 3 tH) /ib with the statistical 
and systematical error, respectively. This value is plotted as the solid point in the lower 
frame of Figure ^. The other data points are from references [58-33,^]. 



The dashed curve shows the energy dependence according to three-body (ppuj) phase 
space and has been normalized to the point directly above this measurement (i.e. ^/s — ^/so = 
0.41 GeV). The solid curve shows the expected behavior when taking into account the finite 
width of the uj meson and the proton-proton final state interaction (FSI) pU| , |B7| . Although 
the energy of this experiment is sufficiently high that pure S wave production is not expected, 
a smooth variation of the cross section with y/s is expected since there are no known baryonic 
resonances with significant branching ratios to pu. Thus, the good agreement of the total 
cross section for u meson production determined here with the curves in the figure, is a 
strong indication that the absolute normalization used here does not introduce an error 
larger than the quoted systematic error. Furthermore, our result is in good agreement with 
the value (45 ± 7) /ib determined in by an interpolation between the existing data. 

The differential cross-section for u meson production has been plotted versus cos Qq^ 
in Figure This distribution has been fit with the sum of the first three even Legendre 
polynomials Pi. The best fit has been obtained with the expression below and shown as the 
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solid curve in the figure. 

^ = (4.0±0.1) •Po + (3.1±0.2) ■P2 + (2.0±0.2) ■P4 (/ub/sr) (5) 

Due to the symmetry of the incoming channel, this angular distribution must have a reflec- 
tion symmetry about cosB^j^^ = 0. Despite the strong variation of the detector acceptance 
with Q^Mi ^his symmetry is observed in the data, providing a very important confirmation 
of the validity of the detector acceptance correction method. The deviations from isotropy 
indicate that partial waves up to L = 2 for the uj meson relative to the pp system are in- 
volved in the production mechanism. Further information such as from polarization degrees 
of freedom are however required in order to make a quantitative measure of the relative 
partial wave amplitudes. 



B. 4> Meson Production 

The distribution is shown in Figure |^ after application of the acceptance correction 

Matrix B (which explicitly includes M^^) and the absolute normalization discussed above. 
This distribution has been analyzed to determine what fraction of the yield is due to non- 
resonant K^K^ production and what fraction has been produced via the resonance. The 
dashed curve shows the estimated non-resonant contribution. The shape of this contribution 
is given by the M^^^ distribution for final states distributed according to four particle phase 
space {ppK^K^). The shape of the resonant contribution is given by the natural line 
shape of the meson folded with a Gaussian to account for the detector resolution. The 
width of this Gaussian is a = 3.3 ± 0.5 MeV, in good agreement with the simulations 
(cr = 3.4 ±0.1 MeV). The total cross section for the pp —>■ ppK^K^ reaction, as well as the 
resonant and non-resonant contributions are summarized in Table |V]. In addition, the total 
cross section for cj) meson production has been included after correcting for the branching 
ratio {Tk+k- l^tot = 0.491) and rounding to two significant digits. 

The pp(f) final state can be defined by two angular momenta, h is the orbital angular 
momentum of the two nucleons relative to each other and I2 is the orbital momentum of the 
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meson relative to the two nucleon CM system. At threshold the pp(j) final state has the 
orbital angular momenta Zi = ^2 = 0. 

The differential cross-section for (p meson production is plotted versus cos in Fig- 
ure]^. This distribution must be symmetric around cosG^^^ = 0. The fact that the observed 
distribution is indeed symmetric about cosG^^ = (within the error bars), as it was al- 
ready observed in Figure |^ for the angular distribution for the u meson, is an additional 
consistency check for the validity of the acceptance corrections. Furthermore, no significant 
deviations from isotropy are seen in these data, indicating that the is predominantly in a 
S wave state relative to the two protons (i.e. /2 = 0). This observation is in good agreement 
with the expectations of Rekalo et al. |^ who claim /2 = up to a meson CM momentum 



of < Mkc which corresponds to an available energy Q ~ 178 MeV (in comparison, this 
measurement is at Q = 83 MeV). 

Moreover, Rekalo suggests that the pp system may be excited to higher partial waves 
at much smaller energies above threshold than needed to excite the meson relative to the 
protons (i.e. /i = l,/2 = should occur at lower energies than li = 0,l2 = !)• To test this 
expectation the proton-proton angular distribution has been evaluated for the events with 
meson production. In the upper frame of Figure ^, the differential cross section is plotted 
as a function of the polar angle under which the protons are emitted relative to the beam, 
measured in the proton-proton reference frame. Here, the observed angular distribution is 
consistent with being isotropic. On the other hand, the proton-proton angular distribution 
exhibits a significant deviation from isotropy, when measured relative to the direction of 
the meson, as presented in the lower frame of Figure ^. This distribution has been 
parameterized with the sum of the first three even Legendre polynomials as listed below and 
shown as the solid curve in the figure. 

^ = (15.0 ±0.9) • Po + (5.1 ±2.2) ■ P2 - (0.7 ±3.6) ■ P4 (nb/sr) (6) 

Odd Legendre polynomials have been omitted due to the refiection symmetry of the final pp 
state about its CM motion direction, which is opposite the momentum in the overall CM 
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frame. These data are well described by using only the lowest two Legendre polynomials, as 
evident by the large uncertainty associated with the P4 coefficient. These results indicate 
that partial waves up to /i = 1 are involved in the proton-proton exit channel. 

The momentum distribution of the particles in the final state is also related to the relative 
partial wave contributions in the ppcf) system. Defining p to be the momentum of a proton 
in the pp reference frame and q to be the CM momentum of the meson, the total cross 
section can be written as the sum of the individual partial wave contributions []69| : 



^ ~ H / Wh,h?dpi^i^. (7) 

Here M/^ /j is the matrix element for a given final state with angular momenta /i,/2; and 
dpi^i^ is an element of the three body phase space which is given by the following formulae. 

dpi.i, - p^'^^V^^'-'dq, (8) 
where the proton momentum is given by 



P=V^L.-^^x^- + ^, (9) 
and the maximum CM momentum of the meson and the available energy are 



^-- = \/^T^' Q = V-s-V^o. (10) 

Assuming that the matrix elements have little variation across the available phase space, 
then the expected differential cross sections as a function of q and p are proportional to 
the variation of the three body phase space with q and p. To illustrate this, the meson 



differential cross section is plotted in the top frame of Figure |T0| as a function of q. The 
dashed curve shows the q dependence of the three body phase space for h = I2 = 0, which 
was normalized to give the smallest relative to the measured data. This curve describes 
the data well. In contrast, the first moment of the dotted curve which denotes the case 
where the meson is in a P wave relative to the nucleons (/i = 0, ^2 = 1) is significantly 
higher than for the data. These results are consistent with the observation from the data 
Figure || that the meson is in a nearly pure S wave state relative to the protons. 
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Although the h = I2 = case describes the data well, a slightly better description of 
the data is generated by including P wave contributions in the proton-proton system (i.e. 
Zi = 0, 1 and /2 = 0). The relative contributions of Zi = and h = 1 has been determined 



by a combined fit (solid curve) to the data in the top frame of Figure [TO together with the 



differential cross section as a function of p, which is shown in the lower frame of Figure IT^. 
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The dashed curve in Figure |l^ denotes the behavior of three body phase space for h = I2 = 0. 
In this case the data appear to have a significantly higher first moment of p, and as a 
result, the best description of the data includes a P wave contribution in the pp system. A 
simultaneous fit to the da /dp and da/dq data (solid curve) yields the following ratio of the 
mean matrix elements for the /i = 1, /2 = to the li = I2 = states: 

0.28 ±0.07. (11) 

Further confirmation that higher partial waves are involved in (f) meson production at 
this beam momentum can be taken from the angular distribution of the daughter kaons from 
(f) meson decay. At threshold (i.e. h = h = 0) the Pauli principle requires the outgoing 
protons to be in a ^5*0 state, and thus the total angular momentum and parity of the system 
must be J'^ = 1^. In this case, angular momentum and parity conservation require the pp 
entrance channel to be in a ^Pi state. Since the orbital projection along the beam direction 
rriL = for the incident plane wave, the angular momentum coupling coefficients require 
the incident pp spin, and hence the outgoing meson spin, to be aligned along the beam 
axis ||7|, p!P| , |55| . Consequently, the angular distribution of the daughter kaons in the meson 



reference frame must display a sin^ distribution relative to the beam direction. At finite 
energies above threshold the spin alignment of the meson is diluted by contributions from 
higher partial waves, thereby modifying the expected angular distribution of the daughter 
kaons. 

The meson spin alignment can be quantified by the spin density matrix. The elements 
of the spin density matrix are related to the emission angles of the kaons from the meson 
decay. After integrating over the azimuthal emission angle (0^") and imposing pu = 
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and pii + Poo + P-i-i = 1-0, the diagonal elements of the spin density matrix are related by 



the following formula to the angular distribution of the daughter kaons (see e.g. |T9|): 



3 



= 2 [P'' ®? + ®?J • (12) 

The differential cross section for the pp ppip reaction has been evaluated and is pre- 
sented in Figure |Tl| as a function of cosG^, based on acceptance Matrix A described in 



Table |. The dotted curve in this figure represents an isotropic distribution and the dashed 
curve corresponds to a sin^ 0^ distribution, which is the expected behavior at threshold due 
to the complete alignment of the meson spin (i.e. poo = 0.0). Within the statistical errors, 
the measured data are not consistent with the sin^ distribution. The solid curve is a fit 
to the data with Eq. |12|, from which the spin density matrix element poo = 0.23 ± 0.04 is de- 
termined. The deviation of poo from the threshold value, together with the result presented 
in Eq. |ll] and the lower frame of Figure ^, indicates a significant admixture of ^-Pi,2 partial 
waves in the outgoing protons for cf) meson production at this beam momentum. {^Pq is 
forbidden in conjunction with I2 = 0, since it would require a l"*" pp entrance channel). 



Although the angular distribution shown in Figure deviates from isotropy, the devia- 
tion is sufficiently small that it justifies as a reasonable approximation the neglect of 0^ in 
acceptance Matrix B, used for the total cross section determination. 

IV. DISCUSSION 

The ratio of the measured total cross sections for the pp —>■ ppcf) and pp —>■ ppuj reactions 
presented in this report is plotted as the filled square in Figure O in comparison to the 



other existing data at higher energies |40-43|. This data point is lower than the other data 



at higher energy. This is primarily due to the different mass of the and uj mesons, and the 
thereby ensuing strong variation of the ratio of available three body phase space volume near 
the (f) meson threshold. These data are compared to a prediction based on a naive application 
of the OZI rule, including the variation of the available phase space (dashed curve). The 
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data point from this measurement is enhanced by roughly one order of magnitude relative 
to the OZI prediction corrected for the available phase space volume. 

The (f)/uj ratio presented here is based on cross sections measured at the same beam 
momentum. In order to reduce uncertainty related to the different phase space volumes, 
partial wave amplitudes, and proton-proton final state interactions, it is useful to compare 
the 0/0; cross section ratio at the same CM energy above threshold. An evaluation based 
on the solid curve in the lower frame of Figure ^ indicates that the cross section for the 
pp — > ppuo reaction is about 8.5 /ib at the same CM energy above threshold as for the cj) 
meson in this measurement (i.e. Q=83 MeV). In this case the (j)/uj ratio would only be 
enhanced by about a factor 5 relative to the OZI rule, in agreement with the higher energy 
data. 



The solid curve in Figure |12| is a calculation from Sibirtsev et al. [jT8| using a one pion 
exchange model and including the proton-proton final state interaction. These calculations, 
which describe the higher energy data well, have an energy dependence similar to the ratio of 
available phase space, and underestimate our point by about a factor of three. Using n — N 
data, Sibirtsev extracts a ratio of the transition amplitudes for ppu to pp(f) production 



R = \M^\/\M^\ = 8.5 ± 1.0. Assuming that R is independent of y/s near threshold, this 
calculation predicts the meson cross section to be (77 ± 16) nb at Q = 83 MeV. A similar 



model by Chung et al. ||70[ including off-shell features of the pion and interference between 
the direct and exchange diagrams predicts about 30 nb at this beam momentum. After 
applying the rather uncertain absolute normalization to our data (see Section P D|) , these 
predictions can be compared to our measured value of 190 ± 14 ± 80 nb. 

Another approach by Nakayama et al. |jT^ explicitly includes not only the mesonic cur- 
rent due to the np ^ (p coupling, but also the nucleonic current where the meson couples 
directly to the nucleon. The observed angular distribution, which is nearly isotropic, indi- 
cates a dominance of the mesonic current in contrast to the cos^ Qf^ distribution expected 
for the nucleonic current. As a result of the dominance of the mesonic current, they can 
not extract a unique value for the coupling constant gNN<t>- In a similar model by Titov et 
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al. |[T^j7T|, our total cross section and differential cross section as a function of Qcm 
also be reproduced without requiring an enhancement of S'^atat over the OZI rule prediction. 
The importance of the correlated vr — p coupling, and the small value of gNN<t> have both 



been predicted by MeiBner et al. |jT5[ to result from large cancellations between intermediate 
kaon and hyperon graphs. 

Although the value oi q^nn is only poorly determined at this point, due to uncertainties 
in the calculations it appears as though our measured data can be explained without invoking 
a large explicit violation of the OZI rule. Nevertheless, the large decay width T^^^^p violates 
the OZI rule itself and requires kaon loop diagrams for a quantitative explanation ||72|. Thus, 
all these solutions require a dominant role of intermediate states with open strangeness, 
thereby indicating that there must indeed be a significant amount of strange sea quarks 
available in reactions involving protons. 

Ellis et al. claim that the enhanced meson production observed in pp reactions proceeds 
dominantly through the "rearrangement" process |[Tl| . To a large degree this is based on the 
observation that the enhanced yield is strongly correlated to the initial spin triplet state 
in pp annihilation |3^. Recent results on np annihilation in flight |73| and pd Pontecorvo 
reactions ||7^ support this hypothesis. According to their model, the "shake-out" process 
should not depend upon the initial spin state, and the "rearrangement" process should 
dominantly occur in the spin triplet state for polarized strange sea quarks, in agreement with 
the data. Following this argumentation, meson production in proton-proton reactions is 
also expected to be strongly correlated with the spin triplet initial state. 

Directly at threshold, the proton-proton entrance channel must be in a ^Pi state due to 
parity and angular momentum conservation, and consequently, the spin of the meson is 
aligned along the beam axis. At Q = 83 MeV we observe the spin density matrix element 
Poo to have a large deviation from the threshold prediction. This deviation is in qualitative 
agreement with the dilution of the spin alignment expected due to the observed contribution 
of the li = 1 partial wave in the exit channel, and suggests that a significant fraction of 
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the meson production at this beam momentum proceeds via the spin singlet ^Sq and 
entrance channels. 

V. SUMMARY AND CONCLUSIONS 

In this paper total and differential cross section values for the production of (p and cu 
mesons in proton- proton reactions at 3.67 GeV/c are presented. The total cross section ratio 
for these mesons is observed to be about an order of magnitude larger than expected from 
predictions based on a naive application of the OZI rule. This enhancement is slightly larger 
than the data measured at higher beam momenta, however significant uncertainty remains 
regarding the relative contributions of different partial waves to the 4> versus lu production 
processes. 

Apparent violations of the OZI rule in pp reactions have sometimes been attributed to a 
significant contribution of intrinsic strangeness to the proton's wave function. On the other 
hand, most of the observed enhancement of the meson yield in pp annihilation can be 
explained in terms of rescattering and loop diagrams. However, the large (f) —>■ up coupling 
also requires kaon loops, and thus the intermediate states are dominated by hadrons with 
strange quark content. Therefore, both interpretations involve a significant contribution of 
strange sea quarks to hadronic reactions involving protons. 

Further information to help determine the origin of the strange sea quarks can be taken 
from polarization observables. For instance, based on a dominance of the rearrangement 
process, the polarized intrinsic strange sea quarks prefer a spin triplet initial state for 
meson production. The differential cross sections presented in this report indicate that 
the proton-proton entrance channel is not in a pure ^Pi state at the beam momentum of 
3.67 GeV/c. If the polarized intrinsic strangeness model is correct, then the ratio 
should increase in direct proportion to the fraction of spin triplet in the initial state. Thus, 
it would be very useful to follow the correlation of the ratio to the spin triplet fraction 
as a function of beam momentum closer to threshold where the triplet fraction must rise. 
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Another sensitive test of the intrinsic strangeness model would be to determine the me- 
son production cross section in proton-neutron reactions. For instance, based on the intrinsic 
strangeness model, Ellis et al. predict the cross section ratio to be On-p^np4>l Opp-*-pp4> ~ 0.25 



near threshold [jTT|. In contrast, meson exchange models [0,0 predict anp^np<i,/ o'pp^pp<f, ~ 5 
near threshold. 

Finally, it would also be very important to determine the uj meson yield and partial 
wave contributions at the same excess energy relative to the threshold as presented here 
for the meson. The angular distributions closer to threshold are needed to disentangle 
the mechanisms involved in the uj meson production |7^, and this would eliminate the 
uncertainty related to the relative contributions of the different partial waves in the ppuo and 
ppcp systems. 



VI. APPENDIX: DATA TABLES 

In this section the differential cross sections presented in the figures above are listed in 
Tables |V1|, |V11| , and |V111| . In addition to the statistical errors quoted here, there are the 
systematic errors and the overall normalization uncertainty discussed above. 
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TABLE I. Kinematic variables and the number of bins per variable associated with the accep- 
tance correction matrices for the ppK'^K~ final state (Matrices A and B), and for the ppX final 
states, where X = r/,aJ. 

Kinematic Variable (Mf,^,^)' (Mf„f)' 6^^ ^^M ^g^M ^i^f Qf 
Matrix A {X = K+K') 4 4 1 1 4 1 10 4 

Matrix B {X = K+K') 1 1 10 1 4 40 1 1 

pp^ppX{X = ri,u) 10 10 20 1 4 - - - 



TABLE IL List of systematic error sources and their effect for the various measured particle 
ratios. 



Systematic Error 


(j)/uj 


u/r] 


K+K~/r] 


(f>/v 


Acceptance Correction 


5% 


5% 


5% 


5% 


uj Background 


15% 


15% 






T] Background 




15% 


15% 


15% 


(p Background 


15% 






15% 


Trigger Bias 


10% 




10% 


10% 


Tracking Efficiency 


10% 


10% 


10% 


10% 


Drift in Electronics 


10% 


10% 


10% 


10% 


K'^K~ Identification 


10% 




10% 


10% 


iT'^Tr~ Identification 


10% 


5% 


10% 


10% 


Total 


32% 


27% 


27% 


32% 
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TABLE III. List of systematic biases and their effect for the various measured particle ratios. 


Systematic Bias 




u/r} 




Non-target Events 


1.05 


1.00 


1.05 


Trigger Bias 


0.93 


1.00 


0.93 


/Ti Tr^^i■n^■^■^p^^^■^^\n 


1 07 




1 07 


a; Identification 


0.93 


1.08 




77 Identification 




0.94 


0.94 


Total 


0.97 


1.02 


0.98 


TABLE IV. Ratios of the total 


meson production 


cross sections for various reaction combina- 


tions including the corresponding 


statistical and systematic uncertainties. 


N.B., the values and 


errors corr('S|)<)])ding lo (lie (')/// nw 


<i c;/^' ratios Imvc 1 


leeu mul(.i]")li(Hi In- 10'^ 





Ratio Statistical Error Systematic Error 



oj/rj 0.37 ±0.02 +0.1-0.08 

(p/r) X 10^ 1.42 ±0.1 +0.45 - 0.34 

4)/lj X 10^ 3.8 ±0.2 +1.2 - 0.9 
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TABLE V. Total production cross section for the reaction pp ppK^K~ at 3.67 GeV/c 
and for the resonant {4> meson) and non-resonant contributions, together with the statistical and 
systematic error, respectively. 

Production CliMnuel Cross Soctiou [//!)] 

Non-resonant K+K' 0.11 ± 0.009 ± 0.046 

(j) K+K- 0.09 ± 0.007 ± 0.04 

Tot al K+R- 0.20 ± 0.011 ± 0.08 

{(j) K+K-) X Ftot/^K+K- 0.19 ± 0.014 ± 0.08 
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TABLE VI. Differential cross sections for uj (left two columns) and (f) (right two columns) meson 
production as functions of cos(9^J^^) and cos(6^J^^), respectively. 



'(©cm) 


da/dO, [//b/sr] 


cos(eSM) 


da/dO, [nb/sr 


-0.9 




-0.9 


8.6 ± 29.5 


-0.7 


3.5 ± 2.0 


-0.7 


15.3 ± 4.5 


-0.5 


3.3 ± 0.4 


-0.5 


17.4 ± 3.3 


-0.3 


2.9 ± 0.34 


-0.3 


16.6 ± 2.0 


-0.1 


3.0 ± 0.29 


-0.1 


17.2 ± 2.0 


0.1 


3.1 ± 0.23 


0.1 


13.3 ± 1.5 


0.3 


2.9 ± 0.13 


0.3 


16.6 ± 1.4 


0.5 


3.3 ± 0.07 


0.5 


14.5 ± 1.2 


0.7 


3.9 ± 0.09 


0.7 


15.4 ± 1.2 


0.85 


5.4 ± 0.1 


0.9 


16.2 ± 1.2 


0.95 


8.0 ± 0.1 
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TABLE VII. Differential cross sections for the pp ppcp reaction as a function of the angle of 
a proton measured in the pp reference frame relative to the beam direction (left two columns) and 
relative to the direction of the (j) nieson (center two columns). The right two columns tabulate the 
differential cross section as a function of the polar angle of the daughter meson, measured in 
the (j) meson reference frame, relative to the beam direction. 



cos(ePp) 


da/dVL [nb/sr] 


cos(*Pp) 


da/d^ [nb/sr] 


cos(ef ) 


da/dQ [nb/sr] 


-0.875 


15.4 ± 2.5 


-0.875 


18.7 ± 2.7 


-0.875 


10.3 ± 2.4 


-0.625 


16.2 ± 1.7 


-0.625 


17.1 ± 3.8 


-0.625 


15.5 ± 1.6 


-0.375 


15.3 ± 1.9 


-0.375 


14.2 ± 2.8 


-0.375 


16.4 ± 1.7 


-0.125 


14.8 ± 1.8 


-0.125 


11.8 ± 2.4 


-0.125 


17.4 ± 1.3 


0.125 


15.9 ± 1.9 


0.125 


12.2 ± 1.9 


0.125 


17.6 ± 1.6 


0.375 


14.1 ± 1.6 


0.375 


14.4 ± 2.2 


0.375 


15.0 ± 2.0 


0.625 


15.8 ± 1.8 


0.625 


14.5 ± 2.4 


0.625 


16.2 ± 1.8 


0.875 


13.5 ± 1.5 


0.875 


18.0 ± 2.9 


0.875 


12.6 ± 2.7 
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TABLE VIII. Differential cross sections for the pp ppcj) reaction as a function of the CM 
momentum of the meson (g, left two columns) and as a function of the proton momentum in the 
proton-proton reference frame (p, right two columns). 



q [GeV/c] da/dq [/xb/(GeV/c)] p [GeV/c] da/dp [/xb/(GeV/c)] 



0.017 


0.011 ± 0.012 


0.0145 


0.0 ± 0.55 


0.051 


0.064 ± 0.015 


0.0435 


0.079 ± 0.040 


0.085 


0.180 ± 0.030 


0.0725 


0.178 ± 0.037 


0.119 


0.379 ± 0.038 


0.1015 


0.199 ± 0.060 


0.153 


0.573 lb 0.052 


0.1305 


0.457 ± 0.058 


0.187 


0.680 ± 0.078 


0.1595 


0.819 ± 0.087 


0.221 


0.923 ± 0.097 


0.1885 


0.961 ± 0.083 


0.255 


1.09 ± 0.087 


0.2175 


1.25 ± 0.12 


0.289 


0.955 ± 0.11 


0.2465 


1.30 ± 0.11 


0.323 


0.585 ± 0.083 


0.2755 


0.926 ± 0.087 


0.357 


0.148 ± 0.098 


0.3045 


0.382 ± 0.065 
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FIG. 1. Schematic layout of the DISTO experimental apparatus, viewed from above. The 
large shaded area represents the effective field region. 
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500 500 

Cherenkov Light [arb. units] 

FIG. 2. Plots of the Cerenkov light output versus particle momentum for negatively(a and b) 
and positively(c and d) charged particles. The figures on the left are inclusive distributions and the 
figures on the right are the same distributions after requiring that the other Cerenkov amplitude 
is consistent with kaon identification and that the four observed particle momenta are consistent 
with 4-momentum conservation for the event hypothesis ppK~^K~ . 
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FIG. 3. Distribution of {M^^f - (A^i^J^ fo^. ^.j^g ^^^^^ hypothesis pp ppK+R-. The 
dashed histogram is an estimate of the background by scahng the inclusive distribution without 
applying kaon Cerenkov requirements. 



41 




(uZf [GeW] 



FIG. 4. Distribution of (M^^^^ )^ versus (M^^^^)^ (upper frame). The lower frame shows 
prominent peaks from the tj and to mesons in the projection onto the (M^-^g)'^ axis after requiring 
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Vs - Vsq [GeV] 

FIG. 5. Exclusive cross section as a function of the total CM energy above threshold for the 
reaction pp — > pprj (upper frame) and of the reaction pp — > ppu (lower frame). The data points are 
referenced in the text. The solid and dashed curves in the upper frame are parameterizations of the 
data in order to estimate the rj meson production cross section at the energy in this experiment, 
which is marked by the vertical dotted line. The filled square in the lower frame is the result of 
this experiment. In the lower frame the dashed curve shows the CM energy dependence of the 
three-body (ppco) phase space volume. The solid curve additionally includes the finite width of the 
CO meson and the proton-proton final state interaction (FSI). 
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-1 -0.5 0.5 1 

COS(0^^ 

FIG. 6. Differential cross-section (in the CM frame) for the pp ppu) reaction as a function 
of cosB^jy^. The sohd curve is a fit to the data with the sum of the first three even Legendre 
polynomials. 
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FIG. 7. Acceptance corrected distribution. The dashed curve is an estimate of the 

non-resonant contribution, and is based on the dependence of four body phase space (ppKK). 

The sohd curve is the sum of the non-resonant and the d) meson contributions. 
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FIG. 8. Differential cross section for cp meson production as a function of cos O^-m- 
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FIG. 9. (Upper frame) Differential cross section for the pp ppcj) reaction plotted as a function 
of cosQ^p, measured in the pp reference frame relative to the beam direction. (Lower frame) 
Differential cross section for the pp ppcf) reaction plotted as a function of cos ^'^p, measured in 
the pp reference frame relative to the direction of the (p meson. The solid line is a fit to the data 
with the sum of the first three even Legendre polynomials. 
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Momentum [GeV/c] 

FIG. 10. Differential cross section for the pp(j) reaction as a function of the CM momentum q 
of the (j) meson (upper frame), and as a function of the proton momentum p in the pp reference 
frame (lower frame). The dashed curves denote the behavior of the three body phase space when 
the (p meson is in a Ss wave state (i.e. h = h = 0) relative to the protons. The dotted curves 
correspond to the phase space distribution for /i = 0, I2 = ^ (upper frame) and for li = 1, I2 = 
(lower frame). The solid lines represent a simultaneous fit to both sets of data presented here with 
the combination of S and P wave contributions in the pp exit channel. 



47 



Q I \ \ \ \ \ \ \ \ \ I \ \ \ \ \ \ \ \ \ I 

-1 -0.5 0.5 1 

COS(0K*) 

FIG. 11. Differential cross section for the pp(j) reaction plotted as a function of the cosine of the 
polar angle of the daughter mesons measured in the (j) meson reference frame, relative to the 
beam direction. The curves represent different assumptions for the (j) meson spin ahgnment along 
the beam axis: zero ahgnment (dotted), full alignment as expected at threshold (dashed), and a 
fit to the data based on partial alignment (solid). 
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FIG. 12. Ratio of the total cross sections for the and ppo; reactions as a function of the 
CM energy above the ^ production threshold. Shown is the value measured in this work (square) 
together with data at higher energies and model calculations described in the text. 
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